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ABSTRACT 

We propose a new method for extracting the non-Gaussian signatures on the isotemperature 
statistics in the cosmic microwave background (CMB) sky, which is induced by the gravitational 
lensing due to the intervening large-scale structure of the universe. To develop the method, 
we focus on a specific statistical property of the intrinsic Gaussian CMB field; a field point in 
the map that has a larger absolute value of the temperature threshold tends to have a larger 
absolute value of the curvature parameter defined by a trace of second derivative matrix of 
the temperature field, while the ellipticity parameter similarly defined is uniformly distributed 
independently of the threshold because of the isotropic nature of the Gaussian field. The 
weak lensing then causes a stronger distortion effect on the isotemperature contours with 
higher threshold and especially induces a coherent distribution of the ellipticity parameter 
correlated with the threshold as a result of the coupling between the CMB curvature parameter 
and the gravitational tidal shear in the observed map. These characteristic patterns can be 
statistically picked up by considering three independent characteristic functions, which are 
obtained from the averages of quadratic combinations of the second derivative fields of CMB 
over isotemperature contours with each threshold. Consequently, we find that the lensing effect 
generates non-Gaussian signatures on those functions that have a distinct functional dependence 
of the threshold. We test the method using numerical simulations of CMB maps and show that 
the lensing signals can be measured definitely, provided that we use CMB data with sufficiently 
low noise and high angular resolution. 
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Introduction 



Determination of the power spectrum of dark matter fluctuations in the observed hierarchical large-scale 
structures of the universe remains perhaps the compelling problem in cosmology. Weak gravitational 
lensing due to the large-scale structure is recognized as a powerful probe of solving this problem as well 



as of constraining the cosmological parameters (Gunn 1967; Blandford et al. 1991; Miralda-Escude 1991 



Kaiser 1992) ) , because it can fully avoid uncertainties associated with the biasing problem. Recently, several 



independent groups have reported significant detections of coherent gravitational distortions of distant 
galactic images (van Waerbeke et al. 200C; Wittman ct al. 2000| ; Bacon, Rcfregier & Ellis 200C; Kaiser] 



Wilson & Luppino 2000; Maoli et al. 2000). On the other hand, the temperature anisotropics in the cosmic 



microwave background (CMB) can be the most powerful probe of our universe, especially of fundamental 
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cosmological parameters (e.g., flu, Sugiyama & Silk 1997). The weak lensing similarly induces distortions 
in the pattern of the CMB anisotropics, and the lensing signatures will provide a wealth of information 
on inhomogeneous matter distribution and evolutionary history of dark matter fluctuations between the 
last scattering surface and present. We then expect that the cosmological implications provided from the 
measurements of lensing effects on the CMB will be very precise, because there is no ambiguity in theoretical 
understanding of the primary CMB physics and about the distance of the source plane. However, it is 
concluded that the weak lensing effects on the CMB angular power spectrum C'i is small (e.g. see |5crjak 



1996 and references therein) , although the detailed CMB analyses need to also take into account the lensing 
contribution. Recently, Seljak & Zaldarriaga (1999) (see also Zaldarriaga fc Seljak 1999| ) developed a new 



method for a direct reconstruction of the projected matter power spectrum from the observed CMB map, 
and showed that it could be successfully achieved if there is no sufficient small scale power of intrinsic CMB 
anisotropies. In this method, the lensing signals can be extracted by averaging quadratic combinations of 
the CMB derivative fields over many independent CMB patches like the analysis to extract the distortion 
effect on distant galactic images, even though the reconstruction maps have a low signal to noise ratio on 
individual patches. 



Excitingly, the high-precision data from the BOOMERanG (de Bernadis et al. 200C; Lange et al 



2000| ) and MAXIMA- 1 ( |Hnany et al. 200C| ; [Balbi et al. 2000| ) have revealed that the measured angular- 
power spectrum Ci is fairly consistent with that predicted by the inflation-motivated adiabatic cold 
dark matter models (also see Tegmark & Zaldarriaga 2000; Hu ct al. 2000| ). The standard inflationary 



scenarios also predict that the primordial fluctuations are homogeneous and isotropic Gaussian (Guth 
fc Pi 1982] ), and then statistical properties of any CMB fields can be completely determined from the 



two-point correlation function C{6) or equivalently C\ based on the Gaussian random theory (Bardeen et 



al. 1986 hereafter BBKS; Bond & Efstathiou 1987 hereafter BE). Taking advantage of this predictability, 



various statistical methods to extract the non-Gaussian signatures induced by the weak lensing have been 
proposed. Bernardeau (1998) found that the lensing alters a specific shape of the probability distribution 
function (PDF) of ellipticity parameter for field point or peak for the Gaussian case as a result of an excess 
of elongated structures in the observed (lensed) map. Although the method could be a powerful probe 
to measure the matter fluctuations around the characteristic curvature scale of CMB, the beam smearing 
effect of a telescope is crucial for the detection because it again tends to circularize the deformed structures. 
Van Waerbeke, Bernardeau & Benabed (2000) then investigated that a statistically correlated alignment 
between the CMB and distant galactic ellipticities could be detected with a higher signal to noise ratio, 
provided that a galaxy survey follow-up can be done on a sufficiently large area. We have quantitatively 
investigated the weak lensing effect on the two-point correlation function of local maxima or minima in the 
CMB map, and it can potentially probe the lensing signatures on large angular scales such as 8 sa 70' that 
corresponds to the matter fluctuations with wavelength modes around A ~ 50/i _1 Mpc (Takada, Komatsu & 



Futamase 200C; Takada & Futamase 2001). Recently, using numerical simulations, it was shown that the 



lensing effect causes a change of normalization factors for three morphological descriptions of the CMB map, 



the so-called Minkowski functionals, against their Gaussian predictions (Bchmalzing, Takada & Futamase 



200C) 



The purpose of this paper is to develop a new simple method for extracting the lensing-induced 
non-Gaussian signatures from the CMB map based on the isotemperature statistics. We then focus on 
specific statistical properties of the intrinsic Gaussian CMB field; a field point that has a larger absolute 
value of the temperature threshold tends to have a larger absolute value of the curvature parameter defined 
by a trace of the second derivative matrix of the CMB field, while the ellipticity parameter similarly defined 
is uniformly distributed independently of the threshold because of the isotropic nature of the Gaussian 
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field. From these features, we expect that the weak lensing causes a larger distortion effect on structures 
of temperature fluctuations around a point with higher threshold. In particular, the lensing can induce 
a coherent distribution of the ellipticity parameter correlated with the threshold owing to the coupling 
between the CMB curvature and the gravitational tidal shear. To extract these characteristic patterns, 
we define three independent functions based on the isotemperature statistics that are obtained from the 
averages of quadratic combinations of the second derivatives of CMB field over isotemperature contours 
with each threshold. As a result, we find the lensing effect on those characteristic functions generates a 
definite functional dependence of the threshold, and therefore the lensing signals could be easily measured 
as a non-Gaussian signature since those functions have very specific shapes in the Gaussian (unlensed) case. 
Using numerical simulations of lensed and unlensed CMB maps including the instrumental effects of beam 
smearing and detector noise, we investigate the feasibility of the method. 

This paper is organized as follows. In §^ we formulate a method for extracting the lensing-induced 
non-Gaussian signatures using the Gaussian random theory for the primary CMB. In we outline the 
procedure of numerical experiments of our method using the simulated CMB maps with and without the 
weak lensing effect. In §[| we present some results in the flat universe with a cosmological constant and 
investigate the detectability of lensing signatures by taking into account the measurement errors associated 
with the cosmic variance and the instrumental effects especially for the future satellite mission Planck 
Surveyo^. In the final section some discussions and conclusions are presented. 



2. Method: Weak Lensing Effect on Isotemperature Statistics 

2.1. Random Gaussian Theory 

In this section, we briefly review a relevant part of the Gaussian random theory developed by BBKS and 
BE for three- and two-dimensional cases, respectively. First, we define the temperature fluctuation field in 
the CMB map as A(0) = [T(8) — Tcmb]/Tcmb- Throughout this paper we employ the flat sky approximation 
developed by BE, and this is a good approximation for our study because the lensing deformation effect 
on the CMB anisotropies is important only on arcminite scales. The Fourier transformation can be 
then expressed as A(0) = J d 2 Z/(27r) 2 A(Z)e^'^, and the statistical properties of the unlensed CMB are 
completely specified by the angular power spectrum Ci defined by (A(i)A(Z')) = (2ir) 2 Ci5 2 (l — I'). 

According to the Gaussian random theory, a certain set of variables Vi(i = 1, 2, N) constructed from 
the CMB field obeys the following joint probability distribution function (PDF); 

P{V) = (27r)^/2|det(M ij )| 

where My is the covariance matrix defined by My = ((uj — (vi))(vj — (vj))), and M _1 and det(My) denote 
the inverse and determinant, respectively. Since we are interested in the lensing distortion effect on the 
isotemperature contours as a function of the temperature threshold, we pay special attention to statistical 
properties of the second derivative field of A, because the local curvature of CMB is probably a good 
indicator of the lensing distortion effect as shown latter. It is then convenient to introduce the following 



(1) 



i 
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variables 

_A A11 + A22 V _A U -A 22 „2A 12 

Z/ = — , A = , Y = , Z = , (J,) 

Co (72 02 CT 2 

where er n is defined by cr 2 = f(ldl/2ir)Cil 2n , Ay = d 2 A/d9 l d8 : > and is the so-called threshold of 
temperature fluctuations. To clarify the physical meanings of X, Y and Z more explicitly, we express them 
in terms of two eigenvalues Ai and A 2 for normalized curvature matrix — Ay /172 as 

X = -(Ai + A 2 ), Y = -2eXcos(2ip), Z = -2eX sin(2y>), (3) 

where e represents the local ellipticity parameter defined by e = (Ai — A 2 )/[2(Ai + A 2 )] and <p denotes the 
relative angle between the principal axis of A !3 - and the 1-axis. By the meaning of equation above, hereafter 
we call X a local curvature parameter around a given field point. Moreover, if the isotemperature contour 
in the neighborhood of local minima or maxima is given by an ellipse of f (01,62) = 0\/b 2 + 2 /a 2 in the 
coordinates of principal axes, the parameter e can be expressed in terms of a and b as e = (a 2 — b 2 ) / [2(a 2 +b 2 )]. 
Hence, Y and Z represent 1- and 2-axis components of the local ellipticity parameter of temperature 
curvature field, respectively. The non-zero second moments of the variables (|^) can be then calculated as 

(v 2 ) = (X 2 ) = 2 (Y 2 ) = 2 (Z 2 ) = 1, (vX) = 7*, (4) 

where 7* = <x 2 /(<7oa2) (although we will also use the same letter 7 for a shear component of lensing 
deformation tensor, we want readers not to confuse 7* and 7) and 7* represents the strength of cross 
correlation between X and v from the relation of 7* = (vX). Equation ([[]) tells us that the joint PDF of 
variables Vi = (y, X, Y, Z) for one field point becomes 

p(u,X,Y,Z)= 2 = exp[-Q], (5) 

(2tt)V 1 -7» 

with 

2Q^ 2 + ( ^~ 7 ;^ )2 +2F 2 + 2Z 2 . (6) 

The important result is that v and X have the non-vanishing cross correlation, and the term of 
exp[— (X — 7* v) 2 1 (2(1 — 7 2 ))] in equation (|5|) physically means that structures around a field point with 
larger absolute threshold tend to have a larger absolute value of the curvature parameter X. In fact, 
this feature is more explicitly clarified by considering the conditional probability distribution for field 
points with a given threshold v. Figure [l] shows the distribution of curvature parameter X subject 
to the constraint that the point has a given threshold u, where the conditional PDF is defined by 
p(X\v) = p(v, X)/p(v) = 1/ yj2n(\ — 7 2 ) exp[— (X — 7*zy) 2 /(2(l — 7 2 ))]. The absence of correlation between 
v and Y or Z is the consequence of the isotropic nature of Gaussian field, more specifically due to the 
isotropic distribution of an orientation angle of ellipticity parameter. 

Using the PDF (||), we define the following three independent functions with respect to temperature 
threshold v% that characterize statistical properties of second derivative fields of CMB along isotemperature 
contours with the threshold i/ t ; 



V x {vt) 


= (<%- 


»t)X 2 ) 


V Y (vt) 


= <*(!/- 




Yziyt) 


= (5(i/- 





V2tt 
1 

2y/2^ 



exp [ ( l_ 7 2 )+7 2 , 2 ] 



ex P 1 — V 



- 5 - 



o.i 



0.4 



0.3 



0.2 - 








-1 







2 



3 



4 



X 



Fig. 1. — The conditional probability distributions of the curvature parameter X for a field point with the 
height v = (solid), v = 1 (dot), v = 2 (short dash), f = 3 (long dash), and v = 4 (dot-dash) are plotted 
for 7* = 0.35. 

where the bracket is defined by (■■•) = JdvdXdYdZ ■ ■ -p(v, X, Y, Z) and can be observationally interpreted 
as average of the considered local quantities performed over all the isocontours in the CMB sky from the 
assumption of large scale statistical homogeneity. All other averages of quadratic combinations of the 
second derivatives such as (XY) = (XZ) = (YZ) vanish because of the isotropic nature of the Gaussian 
field. The functions in equation (0) thus have very specific shapes for the Gaussian case, and we can take 
advantage of this property in order to extract the non-Gaussian signatures on those functions induced by 
the lensing distortion effect. 

2.2. Lensing distortion effect on the isotemperature contours as a non-Gaussian signature 

The CMB photon rays are randomly deflected by the inhomogeneous matter distributions inherent 
in the intervening large-scale structures of the universe during their propagations from the last scattering 
surface to us. Therefore, the observed CMB temperature fluctuation field at a certain angular direction 
6, A(0), is equal to the primary field emitted from the another direction + £(<?) on the last scattering 
surface, A(0 + £), where £(0) is the displacement field. The lensed second derivative field of CMB can be 
then expressed by 



where A is the so-called amplification matrix and is the Kronecker delta symbol. Hereafter, the variables 
with and without tilde symbol denote the lensed and unlensed CMB fields, respectively. The components 




(8) 
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of Aij can be expressed in terms of the local gravitational convergence K and tidal distortion 7 as 

Al3 = { 1 -72 71 l-K + 71 )■ (9) 

In the weak lensing regime the matrix A is always regular, and the variances of k, 71 and 72 are related to 
each other through 

< K 2 )=2< 7l 2 )=2<7 2 2 H^. (10) 

We here have not assumed that k and 7, are Gaussian, and this is a simple consequence of statistical 
isotropy of the displacement field. As shown by the several works using the ray tracing simulations through 
the large-scale structure modeled by N-body simulations, the lensing fields are indeed not Gaussian on 
angular scales of 9 ^ 10' ( Jain, Seljak fc White 2000 ; Hamana, Martel fc Futamase 200C| ), and we will later 



discuss the problem how the non-Gaussian features of k could affect our results. The second moment of the 



convergence, a K , is related to the projected matter power spectrum (Kaiser 1992) 



°l = f^PS) = \h'A f£ fjd X a-\r)W\ XlX ^)P S {k = J- yX ^ , (11) 

where Psik) and P K (l) denote the three-dimensional power spectrum of matter fluctuations and its projected 
power spectrum, respectively, r is a conformal time, X = tq — r, and the subscripts and "r" denote 
values at present and the recombination time, respectively. Hq(= 100/ikm s _1 Mpc" 1 ) and f2 m o denote the 
present-day Hubble constant and energy density parameter of matter, respectively. r(x) is the corresponding 
comoving angular diameter distance, defined as if -1 / 2 sin K 1 / 2 X , X , (— K sinh(— K ) 1 ^ 2 X for K > 0, 
K = 0, K < 0, respectively, where the curvature parameter K is represented as K = (fi m o + £1ao — 1)-Hq 
and fl\o is the present-day vacuum energy density relative to the critical density. The projection operator 
W( X , Xrcc) on the celestial sphere is given by W( X , Xrcc ) — r( Xlcc — X )/r( Xlcc ). As shown later, the effect of 
the finite beam size 9[ w hm of a telescope on a K appears as a cutoff at / >, l sm in the integration of equation 
( pd| ) from the relation of l sin ~ l/#f w hm and thus a K also depends on #f w h m in a general case. Inversely, by 
changing the smoothing scale artificially, we could reconstruct the scale dependence of projected matter 
power spectrum and we will also investigate this possibility. The important result of equation ([ll]) is that 
the magnitude of a K is sensitive to Q m o and particularly to the normalization of matter power spectrum of 
Ps-, which is conventionally expressed in terms of the rms mass fluctuations of a sphere of 8/i _1 Mpc, i.e., 
us ■ Similarly, variances of the second derivative fields of displacement field £i can be calculated as 



(£1,11) = 5 (£,12) - 5 (6,116,2a) = 5 (6,126,22) = y^s 2 (12) 



with 



■s 



2 



Idl 



4 / —l 4 P K (l). (13) 



2tt 



Equation (g) yields the following relations between the lensed (observed) and primary components of 
the second curvature matrix of temperature fluctuations up to the second order of £ : 

X=[(l- k) 2 + 7 2 ] X + 2 7l F + 2 72 Z - (6,ii + 6,22) , 

01 

Y=[(l- k) 2 + 7l 2 - 7I] Y + 2 7i A + ^ (6,11 - 622) , 

Z=[(l- n) 2 - 7l 2 + 72 2 ] Z + 2 72 A + ^6,12- (14) 

cr 2 
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where we have ignored the second order contributions of K^i, 7172 and so on because they vanish after 
the average as a consequence of the statistical isotropy of £j. Note that the weak lensing does not change 
the relations between second moments of these lensed variables compared with the Gaussian cases (^); 
(X 2 ) = 2(Y 2 ) = 2{Z 2 ) fa 1 + l0a 2 K . The equation (0) for Y or Z implies that the lensing effect could 
induce an ellipticity parameter at a certain field point that arises from a coupling between the curvature 
parameter X and the gravitational shear 7 even if the intrinsic ellipticity is zero (Y = Z = 0). Since this 
effect is observable only in a statistical sense, we focus our investigations on the problem how the lensing 
alters statistical properties of the CMB second derivative fields based on the isotemperature statistics. 

Now we present theoretical predictions of lensed functions ([?]) with respect to temperature threshold in 
the observed CMB map. If we assume that the primary CMB fields on the last scattering surface and the 
lensing displacement field due to the large-scale structure are statistically independent, after straightforward 
calculations we can obtain 

\l + 8al){(l-j 2 )+j 2 , 2 }+2a 2 K + ^ i 



V x {vt) = (8{v-v t )X< 



S(v-vt)Y' 
Vy(vt). 



2ir 
1 



ex P 7T 



2V27T 



ex P — 7T 



(1 + 6*1) + 4*1 {(1 - 7*) + 7*V} + 



(15) 



We so far have used the perturbations only for the lensing displacement field £ and thus these equations 
( |l5| ) are valid for an arbitrary threshold v t - Equation ( )l5] ) clearly shows that one of the lensing effects 
on these functions is the change of their normalization factors. The another important effect is that the 
lensing generates a characteristic functional dependence of vt on Vx{vt), VV(^t) an d Vz{vt)- This is as a 
consequence of the coupling between the CMB curvature X and the gravitational tidal shear 7 through the 
intrinsic correlation between i/ t and X, and physically means that the lensing effect distorts more strongly 
the isotemperature contours that have larger absolute threshold. 

In practice it will be rather difficult to discriminate the change of normalization factors caused by 
the lensing in equation ( |l5| ) from the Gaussian case, because measurements of the normalizations in the 
CMB map might be also sensitive to the systematic errors, for example, due to the discrete effect of pixel 
in the map. For this reason, we here focus on the non-Gaussian signatures that have a distinct functional 
dependence of vt , and consider the following observable functions normalized by their values at v t = as a 
deviation from the specific function exp[— v 2 j2\: 



F z (y t ) : 



Vxjvt) 
Vx(0) 

Vyiyt) 
VV(0) 

Vziyt) 

V z {0) 



exp 



ex P — 7T 



exp 



CX P — TT 



(1 + 8^)7^ 



1 - t2 + 10a 2 K 



8^7* 2 



2 J \ + lQ*l-A* 2 Kl r 



exp 



(16) 



where we have neglected the terms including contributions of s 2 in equation ( |i"5j ) because we numerically 
confirmed that the contributions are always small. In the following discussions, these three independent 



functions (16) are compared to the results of numerical experiments. Most importantly, equation ( |16|) shows 
that, although for the Gaussian case in the absence of the lensing we should have Fyfat) — Fziyt) = 
for all vt because of * K = 0, the weak lensing induces distinct non-Gaussian signatures expressed as 
oc v 2 7J exp[— v 2 /2\. Therefore, those two functions can be direct measures of the lensing distortion effect 
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on the isotemperature contours. The property of Fy{v%) = Fz{vt) arises from the statistical random 
orientations of both the CMB ellipticity parameter and the gravitational tidal shear, and we can use 
the relation to distinguish or check the lensing signals against other possible secondary non-Gaussian 
contributions. The lensing effects on Fx, Fy and Fz depend on two parameters of 7* and a K . Then, note 
that 7* is a 'parameter of the primordial CMB anisotropy field, which is not observable, and is related to 
the corresponding direct observable quantity 7, in the lensed CMB map by 7* « 7*(1 — 7er^/2), where 7* 
is defined by 7* = ct 2 /(<7 <72) from a 2 = (A 2 ), a\ — ((VA) 2 ) and erf = ((V 2 A) 2 ). We will therefore have 
to treat 7, as a free parameter in performing the fitting between theoretical predictions (|l6| ) and numerical 
results for the functions Fxiyt), -FV(^t) and Fz(ft) in order to determine a K . We have then confirmed that 
7* is well constrained mainly by Fx(^t)- 



2.3. Effect of filtering 

We so far have ignored the effects of filtering. Actual CMB temperature maps will be observed with a 
finite angular resolution, or the artificial filtering method might be used to reduce the detector noise effect 



( Barreiro et al. 1998 ). The measured temperature map is then given by 

^{e-e s ) = Jd 2 e'w(\e-e'\;0 s )A(e'), (17) 

where W(9; 9 S ) is a window function and throughout this paper we adopt the Gaussian beam approximation 
expressed by W(9\9 S ) = exp[— 9 2 / (29 2 )} / (2tt9 2 ) . For the filtering of a telescope, the smoothing angle 9 S can 
be expressed in terms of its full-width at half-maximum angle #f w hm as 9 S — 9{ w h m /\/8 In 2. The filtered 
lensed temperature field is given by 

A r {0) = Jd 2 0'W{\0 - 0'\;9 S )A (9' + £(0')) . (18) 

Similarly, the filtered second derivatives field of the CMB can be expressed as 

A%(0) = Jd 2 e'W(\0-0'\;9 s )A im (e')A mn (8' + $(0')) A nj (0') 

= Jd 2 0'W(\0 - 0'\;9 S ) [A tij (0') + A, in (0')£n,i(0') + • • •] 

= j-0^W{l-9 s )A l3 {l)e d - e + J *J ^ W(\l + l'\;9 s )A m m^(l')e^+ l '> d + ■ ■ ■ , (19) 

where in the first line of right hand side we have used the part integration and assumed that the surface 
integral is equal to zero for a large CMB survey sky. These equations (|l^) and ( |l9| ) mean that the filtering 
procedure and the lensing effect on the CMB do not commute in a general case. Especially, the last line 
in the right hand side of equation (|l^) shows that the information about a certain mode I' of the lensing 
field £ is coupled to sidebands of the different I modes of the CMB field. The problem of mode coupling 
therefore have to be carefully investigated for accurate measurements of our method. However, since the 
intrinsic CMB anisotropies have a small scale cutoff due to the Silk damping and the directions of the 
CMB curvature and the lensing deformation field are statistically uncorrelated, we could employ a simple 
approximation that the filter function W in the equation ( |l9|) is applied both to the CMB intrinsic field and 
the lensing field independently. The variance of convergence field in equations ( |l6| ) can be then expressed 

by 

- 1 <lH ;W 2 (l;9 s )P K (l), (20) 



(27T) 



2 
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where W(l;9 s ) = exp[— 1 2 0^/2]. Unfortunately, this approximation ( |20| ) might not be so accurate since 
the numerical experiments showed that the magnitude of the convergence field reconstructed by the 
non-Gaussian signatures in the simulated maps is smaller than the value of a K computed by equation (Eo) . 
Therefore, the validity or improvement of this approximation should be further investigated using numerical 
experiments. 



3. Models and Numerical Experiments 

3.1. Cosmological models 

To make quantitative investigations, we consider some specific cosmological models. For this purpose, 
we adopt the current favored flat universe in the adiabatic cold dark matter model, where the background 
cosmological parameters are chosen as f2 m o = 0.3, £1ao = 0.7, h = 0.7, respectively. The flat universe is 
strongly supported by the recent high precision measurements of Ci ( dc Bcrnadis et al. 2000] ; Enany et 



al. 2000). The baryon density is chosen to satisfy f^bo" = 0.019, which is consistent with values obtained 



from the measurements of the primeval deuterium abundance (Buries & Tytler 1998). To compute C/ used 
to make realizations of numerically simulated CMB maps, we used helpful CMBFAST code developed by 
Seljak & Zaldarriaga (1996). As for the matter power spectrum used to compute lensing contributions 
to both numerical and theoretical predictions, we employed the Harrison-ZePdovich spectrum and the 
BBKS transfer function with the shape parameter from Sugiyama (1995). The free parameter is only the 
normalization of the present-day matter power spectrum, i.e., a$. The nonlinear evolution of the power 
spectrum is modeled using the fitting formula given by Peacock & Dodds (1996). 



3.2. Numerical simulations of CMB map with and without lensing 

We perform numerical simulations of the CMB maps with and without the lensing effect following the 
procedure presented in Takada & Futamase (2001) in detail. First, a realistic temperature map on a fixed 
square grid can be generated from a given power spectrum, Ci, based on the Gaussian assumption. Each 
map is initially 60 deg 2 area, with a pixel size of about 0.88 arcmin (= 60 deg/4096). A two-dimensional 
lensing displacement field can be also generated as a realization of a Gaussian process using the power 
spectrum of convergence field, P K , defined by equation (Jj~tj). Note that we have now employed a technical 
simplification that the displacement field is assumed to be Gaussian. As explained, the lensing effect 
is then computed as a mapping between the observed and primordial temperature field expressed by 
A(<?) = A(6 + £). The temperature field on a regular grid in the lensed map is then given by a primordial 
field on a irregular grid using a simple local linear interpolation of the temperature field in the neighbors 
(so-called could- in-cell interpolation). In the case of taking into account the instrumental effects of beam 
smearing and detector noise, we furthermore smooth out the temperature map by convolving the Gaussian 
window function and then add randomly the noise field into each pixel. 

Figure || shows an example of simulated unlensed (left) and lensed (right) CMB maps, where the 
isotemperature contours are also drawn in steps of Ai> t = 0.2. This figure illustrates that the regions around 
crowded contours with higher absolute temperature threshold and larger curvatures are more strongly 
deformed by the lensing. Previous works (Beljak & Zaldarriaga 200C; Zaldarriaga 2000| ) have pointed out 



another but partly similar feature that the power of anisotropies on small scales generated by the lensing is 
correlated with larger scale gradient of the intrinsic CMB field. 
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Fig. 2. — An example of simulated primordial CMB anisotropics map (left) and the lensed map (right). We 
here adopted the adiabatic CDM flat universe model with f2 m o = 0.3, f^o = 0.7, and h = 0.7. The figures 
are on a side of 2 degree and the isotemperature contours are drawn in intervals of Av t — 0.2 (y t = A(0)/oq). 



3.3. The CMB curvature field 

To calculate the second derivative fields of CMB at a certain pixel in the simulated maps, we used a 
method of finite differences between neighboring pixels around the point: 



A >n (i,j) = [A(i - + 1) - 2A(i, j) + A(i + /Sx 2 , 
1 

4<5x 



A Mhj) = — - 1, j - 1) - A(i - 1,3 + 1) - A(t + l,j - 1) + A(i + 1, j + 1)] , (21) 



where 5x is the pixel size and A(i,j) is the local temperature fluctuations at the grid point 



4. Results and Cosmological Implications 

4.1. Numerical results 

In Figure ||, we show the numerical results of the lensed or unlensed three functions ([if]) with respect to 
the temperature threshold v t , which are obtained from each 150 realizations of CMB maps with 60 x 60 deg 2 
area for the filter scale of 3 arcmin. All the curves are plotted at intervals of Av t =0.1. The error bar 
in each bin corresponds to the cosmic variance associated with the measurements and is estimated by 
rescaling the variances obtained from all the realizations by a factor 1/8 when we assume the sky coverage 
of 70% (/sky = 0.7) for a survey of the CMB sky. The figure clearly shows that the lensing deformation 
effect generates a significant functional dependence of v\ on Fy{vt) and Fz(vt) approximately expressed in 
the form proportional to v\ exp[— f t 2 /2]. Especially, the non-Gaussian signatures are pronounced at high 



Fig. 3. — The lensed (left) and unlensed (right) shapes of three functions ([16]) on the isotemperature 
statistics with the filter scale of 3 arcmin and as = 2.0 with respect to the temperature threshold vt, which 
are computed using the numerical simulations of CMB maps. The background cosmological parameters are 
fl m = 0.3, fl\o = 0.7 and h = 0.7. The triangle and round marks in each panel correspond to the averages 
obtained from each 150 independent realizations of the unlensed and lensed maps, respectively. The error 
bars denote the la errors due to the cosmic variance computed for 70% sky coverage of the CMB survey. 
The lensed curves of Fy(v t ) and Fz(vt) clearly show that the lensing generates a significant functional 
dependence of v t approximately expressed by the from of oc v\ exp[— v%/2], whereas the unlensed shapes 
have random errors with both positive and negative values around Fy(v t ) = Fz{ v x) — in each bin of u t . 

absolute threshold such as |z/ t | > 1 as a result of the strong coupling between the gravitational tidal shear 
and the large CMB curvature at such high threshold as explained. For a Gaussian case, Fy{v t ) and Fz{vi) 
should be equal to zero at all bins of v t and thus have both positive and negative large values of the cosmic 
variance in each bin, although the mean values do not exactly converge to zero yet for the number of our 
realizations. These results therefore mean that the non-Gaussian signatures induced by the lensing could 
be significantly distinguished compared to the cosmic variance. In Figure ^ we show the contour map of 
({Fxi — (Fxi))(Fxj — (Fxj)))/(&i&j) calculated from those realizations, where (Fxi) is the mean value and 
<7i the variance of the estimators of Fx{vt) at threshold bin Vi. We have confirmed that the correlations 
for Fy(u t ) and Fz{yt) are similar to the result in this figure. These correlations would be required in 
order properly to quantify the significance of any departure from Gaussian statistics when performing the 
fitting between the numerical results and theoretical predictions. Figure ^ indicates and we have actually 
confirmed that, if wc take the data at intervals of A^t ~ 0.5, the correlation matrix becomes to be very 
close to diagonal in each bin for all the cases we consider in this paper. Taking into account this result, the 
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function Fx (ft) defined by equation (|l6|). Here Fxi and o~i denote the values and variance computed from 
the simulated CMB map in the threshold bin v tl . The contour is stepped in units of 0.2 and the solid and 
dotted contours denote positive and negative values of the coefficients. 



following results are shown in intervals of Af t =0.5. 

Figure ^ demonstrates the results with the filtering scale of 5.5 arcmin computed similarly as in 
Figure ||. One can see that increasing the filtering scale decreases the magnitude of lensing signals for 
Fy(v%) and Fz(vt), because the filtering again tends to circularize the deformed structures in the CMB 
map as pointed out by Bernardeau (1998). However, even in this case of 0f w hm = 5.5' the lensed curves 
of i*V(z/ t ) and Fzivt) remain having the distinct functional dependence of v t compared to the Gaussian 
case of Fy(vt) = Fziyt) — 0. In practice it is important to also take into account the detector noise 
effect on our method. We here assume the instrumental specification of 217GHz channel of the satellite 
mission Planck Surveyor ; the noise level of un — 4.3 x 10~ 6 per a pixel on a side of the FWHM extent 
(0fwhm = 5.5'). The noise field at original fine grid is also convolved with the Gaussian filter of FWHM 
scale to avoid domination of noise spikes at small angular scales ( Barreiro et al. 1998| ) . Figure || shows the 



results. The noise effect reduces the amplitudes of Fx{vt) compared to that in Figure || as a result of the 
change of quantity 7*(= of/ 0002) for the intrisic CMB anisotropics due to the noise. Even in this case, the 
figure clearly shows that the noise level of Planck does not largely change the lensed shapes of Fy and Fz 
compared to the results of Figure ||, although the lensing signals at low thresholds such as \v t \ ^ 1.5 are 
weakened. Importantly, significant non-Gaussian signatures on Fy(i/t) and Fz(vt) still remain compared to 
Fy(vt) = Fz{y%) = 0. This is because our method has so far relied on the normalized observable quantities 
such as v(0) = A(0)/<r and those quantities are more robust against the systematic contributions of the 
detector noise than the CMB fields themselves that are certainly affected by the noise. Figure shows the 
results with as — 2.5 similarly as Figure @. This figure explains that the lensing signals can significantly 
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Fig. 5. — Same as Figure || with a 5.5 arcmin filtering scale and as — 2.0. All the plots are shown in 
intervals of Av t = 0.5 taking into account the result of Figure 0. 



deviate from the unlensed case if the lcnsing effect is adequately large. 



4.2. Cosmological implications from the lensing-induced non-Gaussian signatures 

Another important question we should address is how useful cosmological information on the large-scale 
structure formation can be extracted from the lensing signals onto the isotemperature statistics presented 
in the previous subsection. Since equation (|l6| ) shows that the non-Gaussian signatures depend on 7, and 
a K , we become to consider the problem how accurately the contribution of a K can be determined. As shown 
by Jain & Seljak (1997), the magnitude of a K is sensitive to er§ and Sl m o parameters for the CDM models 



and also depends on the filtering scale 9 S for a realistic case. From equation (20) we find that a K thus 
has following approximate scaling relations for flat universe models around the fiducial model (f2 m o = 0.3, 
SI ao = 0.7 and h = 0.7): 

a K [3'] w 0.157 x (0 fwhm /3')^ a35 x (a 8 /2.0) 1A x (S! m0 /0.3) - 25 , 

cr K [10'] w 9.67 x 10~ 2 x (0 fwhm /lO')" a47 x (as/2.0) 1 - x (r2 mO /0.3) - 10 . (22) 
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Fig. 6. — This figure also includes the detector noise effect for the case shown in Figure ||, where we have 
assumed the noise level of Planck expressed in terms of the variance of noise field per a 5.5 arcmin FWHM 
pixel as a pl ^ — 4.3 x 10~ 6 . Note that the scales of x- and y- axes are same as in Figure |[ 



Furthermore, since more fundamental information is contained in the three-dimensional mass fluctuations, 
we have to take into account the projection effect (Jain & Scljak 1997; Seljak & Zaldarriaga T99S| ; 



Zaldarriaga & Seljak 1999; Takada & Futamasc 2001). The convergence field at 6 £ 10' depends mainly 



on the mass fluctuations with wavelength modes of A ^ 10/i _1 Mpc and the structures distributed in wide 
redshift ranges of z ;$ 10 peaked at z rs 3. The lensing distortion effect on the CMB map can thus be a 
powerful probe of the large-scale structures up to high redshift in principle, which is not attainable by any 
other means. 

Table summarizes the results obtained for the determination of a K with a best-fit and the lcr error, 
which arises from the cosmic variance or also including errors due to the detector noise effect. The 'analytic' 
value of a K in this table is calculated using the approximation for the beam smearing given by equation 
(^0|). We here performed the x 2 -fitting between the numerical results and theoretical predictions for Fx{vt), 
Fy(t't) and Fz{v%)- Note that we have used each data of the functions in the range of —4 < v t < 4 at 
intervals of Ai/ t = 0.5, because the correlation matrix for those data is close to diagonal as explained in 
Figure [I| Then, the theoretical predictions of Fx(vt), Fy{v t ) and Fziyt) are given by two free parameters 
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Fig. 7. — Same as Figure o with erg = 2.5. 



of 7* and er K expressed by equation ( |l6| ) . Most importantly, Table [l] clearly shows that a K estimated from 
the non-Gaussian signatures on Fx, Fy and Fz could be significantly detected with high signal to noise 
ratios compared to the unlensed case. Here, la error for the a K determination corresponds to A\ 2 — 2.3 for 
the x 2 fitting. Figure ^ demonstrates an example of the best-fit results for the noise case with <r 8 = 2.0 and 
$fwhm = 5.5 arcmin. This figure shows that a K is constrained mainly from the numerical data of Fy(vt) an d 
Fz(vt) at | ft | ~ 2 that have relatively small cosmic variances. On the other hand, 7* is well constrained 
only by the data of Fx (ft)- Even if we use the value of 7* directly measured from the lensed simulated 
CMB maps instead of the fitting parameter 7* (see the paragraph under equation (|lq)), it causes only the 
slight change of results in Table |l|. One might then consider a possibility to determine a K by comparing the 
measured 7* from the CMB maps to that of 7* obtained from the fitting of Fx (ft) through the relation of 
7* = 7*(1 — 7(7^/2), but the constraint is much weaker than that of using the non-Gaussian signatures on 
Fy(v%) and Fz(v%). Table [l] also shows that the beam smearing effect is crucial for our method and, in the 
case with #f w h m = 8' and cr 8 = 2.0 the lensing signatures are obscured by the cosmic variance. On the other 
hand, the detector noise effect does not largely affect the results. However, we have to note that the best-fit 
value of a K in all the considered cases is underestimated compared to the analytic value of a K calculated 
by the approximation (^o| ) for beam smearing effect. The possible reasons for the underestimation are 
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Fig. 8. — The results of best-fit between the theoretical predictions and the numerical results for the lensed 
(left: as — 2.0) and unlensed (right) cases with 0f w hm = 5.5 arcmin and the noise effect. The line in each 
panel shows the best-fit theoretical prediction (|16|) (see Table [fl) . 



ascribed to the effect of the discrete pixel in the simulated maps or to the mode coupling between the 



CMB field and lcnsing field caused by the filtering as explained in §2.2. So far it is concluded based on 
the following results that the reason is mainly due to the discrete pixel effect. We have confirmed that the 
best-fit value of 7* even for the unlensed case by our method also underestimates the value of 7, calculated 
by the conventionally used approximation (e.g., Bond fc Efstathiou 198^ ) expressed in terms of Ci and 



#fwhm as 7* = CTi(6>f w hm)/ffo(0fwiim)CTi(0fwiim), wher e of (6 > fwhm ) = / (Idl /2n)l 2l Ci exp[-/ 2 ^] in the context 
of the small angle approximation and 9 S = #f w hm/ Voln2. For example, values of the best-fit and analytical 
prediction for 7* are 0.368 and 0.38, respectively, for unlensed case with 0f w hm = 5.5 and without the 
detector noise effect. For the discrete pixel data of simulated CMB maps, it is generally difficult to perform 
accurate measurements of statistical quantities defined from any derivative fields of CMB compared with 
their analytical predictions. Schmalzing, Takada & Futamase (2000) have also confirmed that it is crucial 
for accurate measurements of Minkowski functionals in a realistic CMB map to take into account the effect 
of discrete pixel, where we have used the interpolation technique. Therefore, we have to further investigate 
the problem how the a K determination from a realistic CMB map performed by our method can reproduce 
its simple analytic prediction, for example, by using the numerical simulations combined to the interpolation 
technique. 
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5. Discussion and Conclusion 

In this paper, we developed a new simple method for extracting the lensing-induced non-Gaussian 
signatures on the isotemperature statistics in the CMB sky and also investigated the feasibility of the 
method using the numerical experiments. Importantly, by focusing on the characteristic three independent 
functions obtained from the averages of quadratic combinations of the second derivatives of CMB field over 
isotemperature contours with each temperature threshold, it was found that the weak lensing generates 
non-Gaussian signatures on those functions that have a distinct functional dependence of the threshold. 
The result is a consequence of the coupling between the gravitational tidal shear and the CMB curvature 
(defined by — [A.n + A ^2]/ ^2) through the intrinsic correlation between the CMB curvature and the 
temperature threshold predicted by the Gaussian theory. By means of the non-Gaussian signatures, it 
can be expected that the lensing signals are extracted from an observed CMB data irrespective of the Ci 
measurements or equivalently the assumptions for the fundamental cosmological parameters. Our numerical 
experiments indeed showed that the method allows us to extract the lensing signals with a high signal 
to noise ratio, provided that we have CMB maps with sufficiently low noise and high angular resolution 



as given by the Planck mission. Recently, Seljak & Zaldarriaga (1999) (see also Zaldarriaga & Seljak 



1999] ) developed a powerful method for a direct reconstruction of the projected power spectrum of matter 



fluctuations from the lensing deformation effect on the CMB maps. The method focuses on the averages 
of quadratic combinations of the gradient fields of CMB over a lot of independent patches in the CMB 
sky like the analysis of the measurements of the coherent gravitational distortion on images of distant 
galaxies due to the large-scale structure. In the method, the lensing signatures on individual patches are 
extracted as differences between the measured statistical measures and their all-sky averages. They showed 
that the reconstruction of input projected matter power spectrum could be successfully achieved if there is 
no sufficient small scale power of intrinsic CMB anisotropics. In this sense, however, the method partly 
depends on the statistical measurements of intrinsic CMB anisotropics, and this is main difference between 
their and our method which we would like to stress. Moreover, our method focuses on the second derivative 
fields of CMB and, therefore, is more sensitive to the amplitudes of the projected matter power spectrum 
on smaller angular scales. Anyway, since the lensing signals in the observed CMB sky are weak, we think 
that some independent statistical methods should be performed to extract them, which could also lead to 
constraints on the projected matter power spectrum at respective, different angular scales. 

Extending our method presented in this paper, one might consider a following possibility to extract the 
lensing distortion effect. BBKS and BE have shown that structures around local higher maxima or lower 
minima of temperature fluctuation field tend to have more peaked shape and be more spherically symmetric 
around the peaks in the Gaussian (unlensed) case. From these features, it can be also expected that the 
weak lensing causes stronger distortion effect on structures around the higher maxima (or lower minima) 
and it might provide us more significant non-Gaussian signatures as a function of the temperature threshold 
of peaks than our method did. However, as quantitatively shown by Bernardeau (1998), the statistical 
measure for the peaks provides the same or lower signal to noise ratio as or than the statistics for a field 
point taking into account the cosmic variance, where he investigated the lensing effect on the probability 
distribution function of ellipticitics around field point or peak. This is because the number density of peaks 
in the CMB sky is not so sufficient for these statistical measurements. For this reason, we prospect a similar 
conclusion for the signal to noise ratio obtained from the measurements of the lensing distortion effect on 
the structures around peaks using our method, although this work will have to be further investigated 
carefully. 



Recently, Schamlzing, Takada & Futamase (2000) have shown that the lensing effect on the Minkowski 
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functionals (the morphological descriptions of the CMB map) appears just as a change of their normalization 
factors against the Gaussian predictions using the numerical experiments. We indeed have confirmed that 
the analytical predictions for the lensed Minkowski functionals done in the similar way as presented in this 
paper give the same conclusion as the numerical results. The result comes from the fact that the lensing 
does not largely change the global topology of the CMB anisotropics in a statistical sense. Likewise, it 
is known that the gravitational potential from which the shear is generated is invariant under the parity 
transformation and the lensing does not induce the so-called 'B-type' polarization defined from combinations 



of the derivative fields of CMB fluctuations ( Seljak fc Zaldarriaga 1999 ; Zaldarriaga fc Seljak 1999 ). These 



results mean that the weak lensing cannot simply generate a new mode of pattern of the CMB anisotropies 
that is absent in the Gaussian case. For these reasons, in this paper we focus on the another information 
on statistical properties of the intrinsic CMB that is useful for the study of lensing effect and can be 
specifically predicted based on the Gaussian random theory. Another issue we should discuss is the possible 
effect on our results caused by non-Gaussian features on the convergence field of the large-scale structure 



at 9 < 10' that are revealed by the ray-tracing simulations (Jain, Seljak & White 2000; Hamana, Martel & 



Futamase 200C), which we have ignored in the numerical experiments of the lensed CMB maps. However, 
since the lensing effect on the CMB can be treated as a mapping, which is expressed as A(9) = A(6 + £), 
and the lensing contributions to the CMB are always coupled to the contributions from the primary CMB 
fields, we prospect that the effect will not change the functional dependence of temperature threshold on 
the lensing- induced non-Gaussian signatures of Fx{vt), -FV(^t) and Fz(y%) expressed by equation (|l6|), even 
if the effect could enhance the magnitudes. 

Undoubtedly, other secondary effects could induce non-Gaussian properties in the observed CMB map. 
The most important effects are the (thermal) Sunyaev-Zel'dovich effect and the foreground contaminations 
such as Galactic foreground and extragalactic point sources. Those effects can be removed to some extent 
by using advantages of their spectral properties, although further reliable investigations should be done for 
any measurements of CMB. Furthermore, in this paper we have presented the lensing-induced non-Gaussian 
signatures on three independent functions, Fx{v t ), Fy^t) and Fz(vt), and therefore we expect that the 
property of Fy(vt) = Fz(vt) can provide us a clue to resolving the lensing contributions from the other 
secondary effects. 



Acknowledgments 

The author would like to thank T. Futamase and J. Schmalzing for fruitful discussions and valuable 
comments. He also acknowledges the useful comments from anonymous referee. He thank U. Seljak and M 
Zaldarriaga for making their CMBFAST code publicly available. He also acknowledges a support from a 
Japan Society for Promotion of Science (JSPS) Research Fellowship. 



REFERENCES 

Bacon, D., Refregier, A., & Ellis, R. 2000, MNRAS, 318, 625 

Bardeen, J. M., Bond., J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 15 (BBKS) 

Balbi, A. et al. 2000, to appear in ApJ Lett., |astro-ph/0005124 

Barreiro, R. B., Sanz, J. L., Martinez-Gonzalez, & Silk, J. 1998, MNRAS, 296, 693 

de Bernadis, P. et al. 2000, Nature, 404, 955 

Bernardeau, F. 1998, A&A, 338, 767 



- 19 - 



Blandford, R. D., Saust, R. D., Brainerd, A. B., & Villumsen, J. V. 1991, MNRAS, 251, 600 
Bond, J. R., & Efstathiou, G. P. 1987, MNRAS, 226, 655 (BE) 
Buries, S. & Tytler, D. 1998 ApJ, 507, 732 
Gunn, J. 1967, ApJ, 147, 61 

Guth, A. & Pi, S.-Y. 1982, Phys. Rev. Lett., 49, 1110 
Kaiser, N. 1992, ApJ, 388, 272 

Kaiser, N., Wilson, G., & Luppino, G. 2000, |astro-ph/0003338 



Lange, A. E. et al. 2000, |astro-ph/0005004 



Hamana, T., Martel, H., & Futamase, T. 2000, ApJ, 529, 56 



Hanany, S. et al. 2000, to appear ApJ Lett., astro-ph/0005125 
Hu, W., Sugiyama, N., & Silk, J. 1997, Nature, 386, 37 
Hu, W. 2000, Nature, 404, 939 



Hu, W., Fukugita, M., Zaldarriaga, M., & Tegmark, M. 2000, |astro-ph/0006463 

Jain, B., & Seljak, U., 1997, ApJ, 484, 560 

Jain, B., Seljak, U., & White, S. D. M. 2000, ApJ, 530, 547 

Maoli, R., V an Waerbeke. L . Mellier. Y., Schneider, P., Jain, B., Bernardeau, F., Erben, T.,& Fort, B. 
2000, |astro-ph/0011251 

Miralda-Escude, J. 1991, ApJ, 380, 1 

Schmalzing, J., Takada, M., & Futamase, T. 2000, ApJ, 544, L83 
Seljak, U. 1996, ApJ, 463, 1 

Seljak, IT., & Zaldarriaga, M. 1996, ApJ, 469, 437 
Seljak, IT., & Zaldarriaga, M., 1999, Phys. Rev. Lett., 82, 2636 
Seljak, IT., & Zaldarriaga, M., 2000, ApJ, 538, 57 
Sugiyama, N. 1995 ApJS, 100, 281 

Takada, M., Komatsu, E., & Futamase, T. 2000, ApJ, 533, L83 

Takada, M. & Futamase, T. 2001, ApJ, 546, 620 

Tegmark, M., & Zaldarriaga, M. 2000, Phys. Rev. Lett., 85, 2240 

van Waerbeke, L., Mellier, Y., Erben, T., Cuillandre, J.-C, Bernardeau, F., Maoli, R., Bertin, E., 

McGracken, H. J., Le Fevre, O., Fort, B., Dantel-Fort Fort, M., Jain, B., & Schneider, P. 2000, A&A, 
358, 30 

van Waerbeke, L., Bernardeau, F., & Benabed, K. 2000, ApJ, 540, 14 

Wittman, D., Tyson, J., Kirkman, D., Dell'Antonio, I., & Bernstein, G. 2000, Nature, 405, 143 
Zaldarriaga, M., 2000, Phys. Rev. D, 62, 063510 
Zaldarriaga, M., & Seljak, U. 1999, Phys. Rev. D, 59, 123507 



This preprint was prepared with the AAS IATjrjX macros v4.0. 



-20- 



Tablc 1. Values of a K from best-fit parameterizations. Errors give the cosmic variance (or also includes 
the instrumental errors caused by detector noise) for 70% sky coverage survey of the CMB map. Here, the 
analytic value of a K is calculated by the approximation (^0|) for the beam smearing effect. 
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